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We respond to the comments expressed by L. D’Alessio in arXiv: 1412.3481 on our work “Floquet 
Fractional Chern Insnlators” [Phys. Rev. Lett. 112, 156801 (2014)]. We confirm the central result 
that the ground state of the effective Hamiltonian is an interacting fractional Chern insulator. 


In this Reply we address the two points of criticism 
regarding our Letter [T] that have been expressed first 
in Ref. [2] and restated in Ref. [3]. These are: i) The 
calculation of the Floquet Hamiltonian of the system 
driven with frequency w misses order terms and ii) 
the assumption that the Floquet bands are filled as in 
time-independent systems is questionable. 

Let us start with the response to statement i). We 
are interested in the effective Hamiltonian governing the 
time evolution of a periodically driven honeycomb lat¬ 
tice system of the form H{t) = Hhop{t) + Lfint- Here 
Hhopit) includes the single particle Hamiltonian and the 
periodic driving and Hint is the interacting Hamiltonian, 
both defined in the first equation of [1] . The stroboscopic 
evolution operator that governs the dynamics of a given 
initial state is, without loss of generality |3] 

U{T + T,t) = eK(r)^-iH„,,T^-ik(T)^ 

where T is the period of the driving. It is written in 
terms of a time periodic kick operator K{t) = K{t + T), 
that encodes all dependence on initial conditions through 
T that averages to zero over one period, and defines the 
T-independent effective Hamiltonian iLeff- In the limit 
where the driving frequency oj = 2TT jT is large compared 
to the other characteristic energy scales of the problem, 
the effective Hamiltonian takes the form (using notation 
of Refs. [HIS]) 

Hes = iJoc + H 1 U .+0 , (2a) 

with Hqui = Hq, (2b) 

00 

= -y^-[Hn,H_n], (2c) 

w n 

n—1 

up to first order in uj~^, where the components Hn 
define the Fourier transform H{t) = J2n=-oo 
Two essential properties of iLeff to order are rele¬ 
vant for our discussion: a) It only includes 
commutators, b) Since the interaction term Hint is 
time-independent in our problem, to order u!~^ the 
only term in which the interaction enters is Hquj = Hq. 
Therefore, there is no interaction correction up to order 


w ^ to the effective Hamiltonian iJeS generically |^. 

Instead of separating the stroboscopic time evolution 
operator as done in Eq. ([^, we could have chosen to 
directly define an evolution operator Hp{t), the so-called 
Floquet Hamiltonian, by writing 

U{T + T,T) = e-^^^^^'>^, (3) 

at the expense of an explicit dependence of Hp[t) on the 
initial time t . The Floquet Hamiltonian Hp{t) reads to 
order uj~^ 

Hpir) = Ho^ + HI,, + O (4) 

with 

Hqui = Hq, (5a) 

n—1 

+e-““"[iL_„,iLo] - e*"“"[H+„,Ho]), (5b) 

In contrast to HeS, we note that Hp (r) includes to order 
u!~^ a) terms that depend explicitly on the initial time 
T and b) interaction corrections of the form discussed in 
Ref. [3] through the commutators [H±n,Ho]. Both a) 
and b) terms generically break point group symmetries 
of the lattice. Importantly, Hp{t) and HeS are unitarily 
equivalent (if all orders of lo are included) 

Hp{t) = ( 6 ) 

As a consequence, the spectra of both Hamiltonians are 
identical. 

The universal time-evolution of the system, indepen¬ 
dent of initial conditions, is thus determined by iLeff- The 
central question that was addressed in our Letter Ref. [T] 
is, whether this time-evolution operator can support a 
topologically ordered Floquet fractional Chern insulator 
(FFCI) as a steady state of the system. The difference be¬ 
tween the effective Hamiltonian used to obtain the results 
in Ref. [1] and HeS are that those terms of [H±n, Hq] that 
add to the non-interacting Hamiltonian were erroneously 
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FIG. 1. Eigenvalue spectrum as obtained with exact diago- 
nalization of Heff in § for a Lx X Ly = 4 X 6 lattice. The 
parameters are the same as in Ref. [I], Fig 3 except that we 
choose a lower amplitude {Ax = Ay = 1.0), rending our pro¬ 
posal even more viable. 
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FIG. 2. a) Flux insertion and b) tower of states as explained 
in Ref. [I] evidencing the fractional and ferromagnetic nature 
of the ground state, respectively. Note that the spectrum in 
a) is symmetric around = n which is not the case when 
H±n] commutators are included. The inset in a) 
shows how the three lowest states permute under flux inser¬ 
tion. 


with spontaneous breaking of spin-rotation symmetry is 
unaffected by this change (see Fig. [^and [^. In fact, the 
gap above the FFCI ground states in a given spin sector 
is found to be larger than in Ref. [T] and the FFCI phase 
is stable at lower amplitudes of the driving field. 

This numerical result does not answer (and did not 
attempt to answer) the question of how to initialize the 
system through the kick operator K (r) in order to obtain 
this desired steady state . It is the power of the decom¬ 
position Q that allows to separate the two problems. In 
contrast, were we to use Hf{t) as Ref. [3] advocates, the 
effect of the kick operators and the time evolution are 
entangled, and thus less transparent to work with. 

This leads us to statement ii) from Ref. [3]: How to 
reach and stabilize the desired FFCI steady state. The 
answer to this question lies in a deliberate control over 
the initialization of the system as well as its coupling 
to a heat bath and goes beyond the scope of Ref. [T]. 
For example, if the bath density of states is suppressed 
at energies that correspond to unwanted transitions, a 
thermal-like population of Floquet could be possible. [51- 
m Once reached, the transitions out of the FFCI state 
may be suppressed for exponentially long times in a limit 
where the driving frequency is large compared to all char¬ 
acteristic energy scales of the undriven system, even in 
absence of a bath. m 

Given the above discussion, we conclude that a) there 
is no correction to order uj~^ in the interaction when cal¬ 
culating HeS, b) this Hamiltonian hosts an FFCI state 
and c) reaching this steady state requires deliberate con¬ 
trol over the initialization of the system and its coupling 
to a bath. In summary, while Refs. [2] and [3] raised im¬ 
portant points about the derivation of effective Hamilto¬ 
nians in time-periodic systems and their relation to the 
actual steady states, we can confirm that our main result 
is unaffected. 
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included in Ref. [T], while they should not be contained 
in a consistent approximation to the effective Hamilto¬ 
nian to order l/oj. However, the presence or absence of 
these terms does not affect qualitatively the main result 
of Ref. [T], namely an affirmative answer to the above 
question. It was pointed out in Ref. [7] that, even though 
the terms [FI±n,i7o] appear at order l/oj in the Floquet- 
Magnus expansion, they affect the spectrum only at order 
l/w^. As a numerical confirmation, we have exactly di¬ 
agonalized as given by Eq. (H) with the terms up 
to order i.e. without the r-dependent terms. We 
found that the nature of the ground state being an FFCI 
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